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The  study  of  start-up  performance  for  a  direct  internal  reforming  molten  carbonate  fuel  cell  (DIR-MCFC) 
system  is  presented.  Since  a  kW-class  stack  is  assembled  with  an  additional  preheating  design,  the 
improvement  of  start-up  behavior  is  conducted  to  find  the  proper  operating  strategy.  For  a  cold  start¬ 
up  fuel  cell  system,  both  start-up  delay  and  inverse  response  are  strictly  detected.  When  the  optimum 
operating  strategy  is  determined  by  solving  the  steady-state  optimization  algorithm  subject  to  stack  tem¬ 
perature  constraint,  the  rapid  system  start-up  as  well  as  the  maximum  power  output  can  be  achieved 
simultaneously. 
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1.  Introduction 

The  molten  carbonate  fuel  cell  (MCFC)  systems  are  being  con¬ 
sidered  as  one  of  promising  solutions  in  new  energy  technologies, 
because  it  could  provide  the  higher  cell  voltage,  the  better  ther¬ 
mal  efficiency  and  the  longer  life  than  other  fuel  cell  systems  [1]. 
Hydrogen  is  ideal  fuel  for  a  general  class  of  MCFC  systems,  but  the 
hydrogen  technology  in  storage,  transportation  and  economic  effi¬ 
ciency  are  not  mature  yet  in  the  commercial  market.  At  the  present 
stage,  alternative  fuel  sources  instead  of  pure  hydrogen  through  the 
hydrocarbon  steam-reforming  process  were  effective  approaches 
[2].  For  high  temperature  MCFC  systems,  methane  is  often  treated 
as  the  primary  fuel.  Through  internal  reforming  [3]  or  pre-reformer 
[4],  the  MCFC  system  could  provide  the  high  energy  efficiency  and 
low  pollutant  emissions  compared  to  traditional  fossil  fuel-based 
energy  systems  [5].  Recently,  Sundmacher  et  al.  [6]  introduced  sev¬ 
eral  approaches  for  an  industrial-scale  MCFC  plant  by  virtue  of 
modeling,  analysis,  optimization  and  control  techniques.  In  gen¬ 
eral,  this  industrial-scale  system  start-up  usually  took  hours  and 
days  due  to  the  requirement  of  too  high  operating  temperature 
and  specified  heating  demands. 

As  for  determining  the  proper  operating  manner  and  examining 
quantitative  response  of  the  system  to  different  inlet  perturbations, 
the  mathematical  modeling  and  simulation  for  a  class  of  MCFC 
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systems  have  been  investigated  [7-9].  Obviously,  the  steady-state 
and  dynamic  characteristics  of  two-dimensional  model  through 
numerical  techniques  and  model  reduction  method  are  used  to 
establish  a  simple  simulation  model  according  to  various  operating 
conditions.  Based  on  the  mathematical  modeling,  the  temperature 
control  of  a  specified  MCFC  model  for  performance  improvement 
was  introduced  [10],  the  optimum  operating  policy  was  deter¬ 
mined  by  solving  a  dynamic  optimization  problem  with  respect  to 
system  efficiencies  [11],  and  the  evaluation  of  energy  and  exergy 
efficiencies  of  an  integrated  power  generation  has  been  addressed 
in  [12]. 

To  study  the  start-up  operation  of  high  temperature  fuel  cell 
systems,  the  dynamic  modeling  and  analysis  of  different  PEM  fuel 
cell  systems  were  widely  studied  due  to  the  prospective  automotive 
market  [13,14].  For  most  solid  oxide  fuel  cell  (SOFC)  stack  systems, 
the  typical  start-up  procedure  would  depend  on  combustor  and 
heat  exchanger  for  heating  of  the  air  and  stack.  The  modeling  and 
thermal  dynamics  of  the  SOFC  stack  during  start-up  process  have 
been  investigated  in  [15].  Zhou  et  al.  [16]  showed  that  the  start¬ 
up  performance  of  MCFC  stack  depended  on  different  temperature 
and  pressure  effects.  However,  the  dynamic  modeling  of  the  MCFC 
start-up  process  according  to  various  operating  manner  is  seldom 
investigated. 

Since  the  fast  start-up  design  on  high  temperature  fuel  cell 
systems  would  face  some  challenges  due  to  different  preheating 
manners  and  operating  conditions,  the  dynamic  modeling  of  a 
jacket-integrated  DIR-MCFC  stack  is  studied  initially.  According  to  a 
simple  heat  exchanger  structure  as  a  preheater  design,  the  dynamic 
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Nomenclature 

^cell 

effective  cell  area,  cm2 

cs 

P 

average  stack  heat  capacity,  J  g-1  K-1 

rhj 

heat  capacity  of  jacket,  Jg-1  K_1 

F 

Faraday  constant  (=96485  C  mol-1 ) 

hf1 

inlet  specific  enthalpy  of  species  i  into  the  anode, 

Jmol-1 

hf 

inlet  specific  enthalpy  of  species  j  into  the  cathode, 

Jmol-1 

lstack 

stack  current  density,  mA  cm-2 

ncell 

number  of  cells  in  the  stack 

NJ1,  N[cn 

inlet  molar  flow  rate  at  the  anode/cathode,  mol  s-1 

Na,Nc 

total  moles  at  the  anode/cathode,  mol 

Pa,i 

partial  pressure  of  species  i  at  the  anode,  Pa 

Pcj 

partial  pressure  of  species  j  at  the  cathode,  Pa 

R 

universal  gas  constant  (=8.314J  mol-1  K-1 ) 

ms 

total  stack  mass,  g 

mhj 

mass  of  steam  in  the  jacket,  g 

™hj 

mass  flow  rate  of  steam  in  the  jacket,  g  s-1 

Rt 

thermal  resistance  of  the  stack,  KW-1 

^amb 

ambient  temperature,  K 

UA 

product  of  the  overall  heat  transfer  coefficient  and 

heat  transfer  area,  W  K-1 

vin 

*a,i 

inlet  mole  fraction  of  species  i  into  the  anode 

vin 

*c,j 

inlet  mole  fraction  of  species  j  into  the  cathode 

Heating  Jacket 

|  Steam 


Fig.  1.  The  run  schematic  of  the  DIR-MCFC  stack  start-up  operation. 


simulation  for  the  DIR-MCFC  stack  would  be  conducted  to  study  the 
system  start-up  operation.  For  a  cold  start-up  system,  both  start¬ 
up  delay  and  inverse  response  would  strongly  increase  difficulties 
on  the  development  of  control  system.  To  address  the  rapid  sys- 


time(s) 

Fig.  2.  No  preheating:  the  stack  temperature  responses  under  different  steam-to- 
carbon  ratios. 


tern  start-up  as  well  as  the  maximum  power  output,  the  optimal 
operating  condition  depends  on  undetermined  jacket  temperature 
and  steam-to-carbon  ratio.  Through  the  steady-state  optimization 
algorithm  subject  to  the  upper  limit  of  stack  operating  tempera¬ 
ture,  the  optimum  air  flow  rate  is  obtained  which  is  verified  by 
simulation  results. 


2.  System  description 

Inspired  by  a  250  kW  class  of  MCFC  stack  where  the  anode  and 
cathode  gas  recirculation  manners  were  utilized  to  recover  waste 
heat  as  well  as  improve  the  power  efficiency  [17],  the  new  heat 
integration  design  for  the  MCFC  stack  system  is  depicted  in  Fig.  1. 
At  the  room  temperature  (298  K),  the  cold  mixture  of  methane  and 
water  is  preheated  by  a  traditional  heat  exchanger  (FIEX)  and  then 
the  heated  mixture  with  specified  steam-to-carbon  ratio  is  directly 
fed  into  the  internal  reforming  unit.  In  the  meanwhile,  the  stack 
temperature  is  perheated  by  a  jacket  with  inlet  hot  steam  until  the 
start-up  temperature  for  the  reforming  process  is  achieved  at  673  K 
at  least.  Afterwards,  the  feed  would  be  converted  into  hydrogen  rich 
gas  mixture  via  the  following  two  reactions: 


Fig.  3.  For  a  hot  start-up  system  under  different  air  flow  rates:  (a)  stack  temperature  responses  and  (b)  corresponding  stack  voltage. 
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Table  1 

Kinetic  parameters  of  reforming  reactions. 


Ar  =  3.71  X  1017 

AEr  =  240.1  kj  mol”1 

7W  =  5.43 

A£w  =  67.13kJmol-1 

BCh4  =  6.65  x  10“9 

AHCh4  =  -38.28  kj  mol-1 

Bh2o  =  1.77x105 

AHh2o  =  88.68  kj  mol”1 

BHi  —  6.12  x  10“14 

AHC0  -70.65  kj  mol”1 

Bco  —  8-23  x  10“10 

AHh2  =  -82.9  kj  mol-1 

Fig.  4.  For  a  cold  start-up  system:  the  responses  of  stack,  jacket  temperatures  and 
stack  voltage. 

(i)  Steam/methane  reforming  reaction: 

CH4  +  H2oA+CO  +  3H2,  AH°98  =  206  kj  mol-1  (1) 

(ii)  Water-gas  shift  reaction: 

CO  +  H2oA^C02  +H2,  AH^gg  = -41  kjmol-1  (2) 

where  both  rate  equations  rR  and  rw  are  shown  in  Eqs.  (Al )  and 
(A2).  Under  an  external  current  demand,  the  half-cell  reactions  take 
place  in  anode  and  cathode  electrodes,  respectively  (Table  1), 


Anode: 

H2  +  C032“  -*  C02  +  H20  +  2e~ 

(3) 

Cathode: 

lo2  +  C02  +  2e-^C032- 

(4) 

The  overall  cell  reaction  occurs  in  a  homogenous  layer  to  pro¬ 
duce  electricity  and  heat  simultaneously.  The  anode  exhaust  gas 
is  fully  oxidized  with  the  excess  air  in  a  catalytic  combustor  and 
then  the  produced  carbon  dioxide  with  oxygen  is  fed  into  the  cath¬ 
ode  channel.  The  outlet  jacket  steam  flows  into  the  external  HEX  to 
warm-up  the  fresh  mixture  of  methane  and  water.  In  our  proposed 


time(s) 


time(s) 


Fig.  5.  For  a  cold  start-up  system  under  different  steam-to-carbon  ratio:  (a)  stack 

temperature  responses  and  (b)  corresponding  stack  voltage. 

system  configuration,  some  general  assumptions  and  start-up  con¬ 
ditions  are  noticed  as  follows: 

(i)  Gas  mixtures  at  inlet/outlet  flows  of  the  stack  are  ideal  mix¬ 
tures. 

(ii)  Thermodynamic  properties  follow  the  ideal  gas  law. 

(iii)  Energy  accumulates  only  in  the  stack  solid  mass. 

(iv)  The  outlet  flow  temperatures  of  the  DIR-MCFC  stack  are  equal 
to  the  stack  temperature,  and  the  outlet  jacket  temperature  is 
equal  to  the  jacket  temperature. 

(v)  There  is  no  phase  transition  between  water  and  steam  in  the 
inside  of  DIR-MCFC  stack. 

(vi)  There  are  no  concentration,  velocities  and  temperature  gradi¬ 
ents  perpendicular  to  the  direction  of  gas  flow. 

(vii)  The  anode  and  cathode  are  well  mixed  volumes  with  inter¬ 
change  of  carbonate  ions  through  the  electrolyte  matrix 
separating  the  two  sides. 

(viii)  The  hot  steam  is  fed  into  the  heating  jacket  and  the  exhaust 
steam  flows  into  the  HEX. 

(ix)  During  the  stack  start-up  period,  the  heating  jacket  instead 
of  HEX  acts  to  heat  the  mixture  of  methane  and  water  at  the 


Table  2 

Heat  capacity  coefficients  for  species  from  298  to  1500  K  at  constant  pressure. 


Species 

href  (kcalmol-1) 

a  (cal  K-1  mol-1 ) 

b  ( x  1 0-3  cal  K_1  mol-1 ) 

c  ( x  1 0-8  cal  K-1  mol-1 ) 

h2 

0.000 

6.947 

-0.200 

0.481 

o2 

0.000 

6.148 

3.102 

-0.923 

H20 

-57.800 

7.256 

2.298 

0.283 

co2 

-94.054 

6.214 

10.396 

-3.545 

CO 

-26.417 

6.420 

1.665 

-0.196 

ch4 

-17.895 

3.381 

18.044 

-4.300 
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Fig.  6.  For  a  cold  start-up  system  under  different  jacket  temperatures:  (a)  stack 
temperature  responses  and  (b)  corresponding  stack  voltage. 


anode,  and  meanwhile  the  air  is  heated  with  the  heating  jacket 
at  the  cathode. 

(x)  For  the  sake  of  simplicity,  the  dynamics  of  external  devices 
including  HEX  and  catalytic  combustor  are  omitted. 

(xi)  The  overall  system  is  a  continuous  process  from  start-up  stage 
to  steady-state  operation. 


According  to  above  statements,  the  dynamic  modeling  of  DIR-MCFC 
is  introduced  in  Appendix  A.  The  state-space  representation  of  the 
DIR-MCFC  system  is  described  by: 


$  =  F(M) 
y  =  m 


(5) 


where  the  state  '£,  =  [xai,  xC  J ,  T°ut,  T]t  represents  species  i  (Table  2) 
at  the  anode,  species  j  at  the  cathode,  jacket  temperature  and  stack 
temperature.  The  non-linear  vector  function  F  is  constructed  from 
Eqs.  (A1)-(A13).  8  =  [x”CH4,  xj.n02,  T”f  represents  the  inlet 

mixture  of  methane  and  water  at  the  anode,  air  (oxygen)  flow  rate 
at  the  cathode  and  inlet  jacket  temperature,  respectively.  The  cell 
voltage  is  denoted  as  the  output  variable  y  =  Vcell,  and  the  corre¬ 
sponding  output  function  h(§)  is  formulated  by  Eqs.  (A14)-(A16). 


Remark  1.  Basically,  the  state-space  model  in  Eq.  (5)  is  used 
to  describe  the  dynamic  behavior  of  jacket-integrated  DIR-MCFC 
stack  regardless  of  the  system  start-up.  If  the  preheating  pro¬ 
cedure  is  neglected,  and  the  temperature  and  pressure  of  the 
experiment-based  DMFC  reliably  range  from  673  to  923  K  and 
1-3  atm,  respectively,  then  the  system  responses  with  respect  to 


Fig.  7.  For  a  cold  start-up  system  with  specified  e  =  2 :  (a)  stack  temperature 
responses  under  different  jacket  temperatures  and  prescribed  air  flow  rate  and  (b) 
corresponding  stack  voltage. 


the  current  demand  are  investigated  by  the  following  simulation 
results 


•  If  the  inlet  flows  of  the  system  is  heated  up  to  the  lower  limit  of 
stack  operating  temperature,  i.e.  Tmin  =  673  K,  Fig.  2  depict  the 
stack  temperature  responses  with  respect  of  different  steam-to- 
carbon  ratios,  £  =  ^anH2o/xancH4  at  a  constant  air  (02)  flow  rate 
NJ.nxJ.no2  =  4mols-1.  Since  the  steam-reforming  of  methane  is 
an  endothermic  process,  the  conversion  of  methane  is  limited 
by  chemical  equilibrium  and  temperature  effect.  Therefore,  the 
large  steam-to-carbon  ratio  induces  the  lower  stack  temperature 
at  steady-state  due  to  a  large  amount  of  water. 

•  In  this  case  of  dynamic  model  of  the  MCFC  stack,  the  desired  oper¬ 
ating  temperature  is  923  K.  At  a  temperature  lower  than  773  K, 
the  cell  performance  drops  significantly,  while  a  temperature 
higher  than  923K  accelerates  material  corrosion  and  power  loss. 
Under  a  given  steam-to-carbon  ratio  e  =  2,  Fig.  3(a)  shows  that  the 
increase  of  air  flow  rate  can  increase  stack  temperature.  Accord¬ 
ing  to  the  upper  limit  of  the  stack  temperature,  i.e.  Tmax  =  923  K, 
Fig.  3(b)  demonstrates  significantly  voltage  drop  in  the  cell  circuit 
while  the  air  flow  rate  is  more  than  4  mol  s-1 . 


Apparently,  the  operating  temperature  strictly  affects  the  stack  per¬ 
formance,  and  the  improper  steam-to-carbon  ratio  or  air  flow  rate 
induces  significant  voltage  drop. 
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Fig.  8.  For  a  cold  start-up  system  with  specified  s  =  3:  (a)  stack  temperature  responses  under  different  jacket  temperatures  and  prescribed  air  flow  rate  and  (b)  corresponding 
stack  voltage. 


3.  System  start-up  performance 

Actually,  the  preheating  procedure  is  essential  when  the  stack 
temperature  is  lower  than  the  normal  operating  temperature.  To 
address  the  start-up  model  of  the  jacket-integrated  DIR-MCFC 
model,  the  state-space  equation  in  Eq.  (5)  is  reduced  to 

f)  =  q(q,  8),  0  <  t  <  tr  (6) 

where  the  warm-up  period  tr  >  0,  the  state  q  =  [T°ut,  T]t  and  the 
vector  function  q  by  Eqs.  (A10)  and  (A12)  is  rewritten  as: 


Remark  2.  According  to  above  start-up  model,  there  are  no 
electrochemical  reactions  in  the  cell  during  a  warm-up  period. 
The  heating  jacket  carries  out  the  preheating  procedure  includ¬ 
ing  heating  mixture  of  methane  and  water  and  air  flow  rate  at  the 
anode/cathode  as  well.  However,  the  warm-up  period  is  uncertain 
due  to  variety  of  inlet  conditions.  In  Eq.  (7),  the  steam-to-carbon 
ratio,  air  flow  rate  and  inlet  jacket  temperature  dominate  the  tran¬ 
sient  behavior  of  the  start-up  model.  Furthermore,  the  preheating 
and  start-up  operation  via  the  dynamic  simulation  are  discussed  as 
follows: 


Q  = 


mhj 


Qheat 


1_  ( yin  _ yout\  _ 

mhj  1  hj  '  mh]CpJ 


NfTntsA  (Jf  -  ftp  +^Ef6ScX5  (5  -  ~hj) 

msCp  +  (Ea,h  +  Ec,h)/T 


hi"  -  hi)  -  {T  -  Tamb)/Rt  +  Qheat 


(7) 


n^ref 


and  its  initials  t/( 0)  = 


W.  Wu  et  al.  /  Journal  of  Power  Sources  195  (2010)  6732-6739 


6737 


Fig.  9.  For  a  cold  start-up  system  with  specified  s  =  5:  (a)  stack  temperature 
responses  under  different  jacket  temperatures  and  prescribed  air  flow  rate  and  (b) 
corresponding  stack  voltage. 


Table  3 

Optimum  air  flow  rates. 

Inlet  condition 

e  =  2,  T“?  =  673  K 
s  =  2,T™  =  723 1< 

8  =  2,  T™  =  723  K 
£  =  3,  T“?  =  673  K 
£  =  3 ,7™  =  723  K 
£  =  3  ,T“  =  723  K 
£  =  5,  T“?  =  673 1< 

£  =  5,  T™  =  723  K 
£  =  5  ,Tjj}  =  723  K 


Air  flow  rate  u( 0)  (mol  s_1 ) 

4.97 

4.71 

4.4 

5.34 

5.07 

4.78 

6.67 

6.41 

6.12 


By  above  simulation  study,  the  stack  start-up  delay  and  the  stack 
voltage  are  all  affected  by  varieties  of  carbon-to-steam  ratio  and 
inlet  jacket  temperature.  However,  a  rapid  system  start-up  opera¬ 
tion  may  induce  a  low  power  output  or  voltage  drop. 


4.  Optimum  operating  conditions 


Referring  previous  simulation  tests,  the  high  stack  temperature 
can  improve  the  power  performance,  but  the  too  high  stack  tem¬ 
perature,  i.e.  T>Tmax,  would  induce  significant  voltage  drop.  We 
found  that  the  variety  of  inlet  conditions  including  carbon-to-steam 
ratio  (£),  inlet  jacket  temperature  (T£)  and  air  flow  rate  (N^xJT^ ) 
strongly  affect  the  stack  temperature  in  a  steady-state  operation  as 
well  as  stack  start-up  delay  in  a  start-up  operation.  To  address  a 
rapid  system  start-up  operation,  i.e.  small  start-up  delay,  and  also 
achieve  the  desired  stack  temperature  in  a  steady-state  operation, 
the  optimum  operating  strategy  is  an  essential  approach.  Moreover, 
a  combination  of  jacket-integrated  DIR-MCFC  and  start-up  model 
is  described  by: 


f  ff  =  q(ij,S ),  0<  t<  tr 

\$  =  F&S),  t>t 


•  When  initial  temperatures  of  inlet  flows  are  298  K,  Fig.  4  shows 
profiles  of  the  stack  temperature,  jacket  temperature  and  stack 
voltage,  respectively,  under  the  inlet  conditions  including  T™  = 

723  K,  6  =  2,  Nlcnx™ q2  =  4mols-1  and  istack  =  100 Am-2.  During 
the  start-up  period,  the  warm-up  period  is  estimated  while  the 
stack  temperature  approaches  to  Tmin  (673  K)  and  the  system 
start-up  delay  is  equivalent  to  tr.  Since  the  jacket  is  a  simple  heat 
transfer  device,  it  is  treated  as  the  preheater  in  the  start-up  period 
but  it  may  become  a  cooling  device  while  the  stack  tempera¬ 
ture  is  higher  than  the  jacket  temperature.  Therefore,  the  inverse 
response  appears  due  to  a  change  of  sign  of  temperature  differ¬ 
ence  between  jacket  and  stack  at  the  start-up  and  steady-state 
stages,  respectively. 

•  Under  T™  =  723 1< and  N^x1^  =  3mols-1,  Fig.  5(a)  shows  the 
responses  of  stack  and  jacket  temperatures  under  various  steam- 
to-carbon  ratios.  The  corresponding  stack  voltage  responses  are 
depicted  in  Fig.  5(b).  Notably,  the  small  value  of  e  can  increase 
system  temperature  and  enhance  the  stack  voltage  at  stead- 
state  operation,  but  the  corresponding  start-up  delay  is  obvious. 
Apparently,  different  carbon-to-steam  ratios  may  encounter  the 
trade-off  between  small  start-up  delay  and  large  power  output. 

•  Under  6  =  2  and  N'J'x1^  =  3mols-1,  Fig.  6(a)  shows  the  stack 
temperature  responses  under  various  inlet  jacket  temperatures. 
Notably,  the  increase  of  inlet  jacket  temperature  can  gradually 
increase  the  stack  temperature  and  meanwhile  reduce  the  start¬ 
up  delay  depicted  by  the  zoom-out  plot  in  Fig.  6(b).  However, 
the  highest  jacket  temperature,  e.g.  T ™  =  723  K,  would  induce 
voltage  drop  at  steady-state  due  to  the  steady-state  stack  tem¬ 
perature  more  than  Tmax  (923 1<). 


where  r\{ 0)  = 


Kj’Tref 


and 


?(fr)  =  Kni,x-,r°ut,Tmin]T. 


Obvi¬ 


ously,  the  inlet  conditions  of  S  would  dominate  the  overall  dynamic 
behavior  and  the  start-up  delay.  Since  a  rapid  start-up  operation 
may  cause  the  degradation  of  stack  performance  due  to  improper 
carbon-to-steam  ratio  and  too  high  stack  temperature,  the  opti¬ 
mum  start-up  operating  policy  is  required  under  a  steady-state 
optimization  framework.  Moreover,  the  minimization  of  the  fol¬ 
lowing  objective  function  subject  to  system  constraints  is  shown 
by: 


min  J  =  (TSs  —  Tmax) 

U(0) 


(9) 


s.t., 

0  =  F($ss,u(0)) 
Tmin  -  Tss  <  Tmax 


(10) 


where  the  initial  condition  u(0)  e  1 6,  N^xf^,  7™  j  and  §Ss  repre¬ 
sents  the  stack  steady-state  value.  Notably,  it  is  typically  non-linear 
programming  (NLP)  problem  to  find  the  solution  u(0).  The  solver 
from  Matlab  Optimization  Toolbox®  is  used. 


Remark  3.  According  to  the  optimization  algorithm,  one  of  inlet 
conditions  can  be  obtained  via  the  numerical  computation  in  the 
Matlab  environment.  The  off-line  iteration  is  used  to  find  the  mini¬ 
mization  of  single  objective  in  a  steady-state  operation,  i.e.  lim  T  = 

t — >oo 

Tmax-  In  Table  3,  optimum  air  flow  rates  u(0),  are  obtained  according 
to  various  inlet  jacket  temperatures  T ™  =  673,  723,  773  K,  and  a 

given  steam-to-carbon  ratio  6  e  {2,  3,  5}.  The  dynamic  responses 
of  stack  temperature  and  stack  voltage  are  depicted  in  Figs.  7-9, 
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time(s) 


time(s) 


time(s) 


Fig.  10.  The  responses  of  molar  flow  rates  of  each  component  in  the  stack  at  different 
operating  conditions:  (a)  (e,  ,  T™)  =  (2,  4.97,  673);  (b)  (e,  ,  T™)  = 

(3,  5.34,  673);  (c)  (e,  ,  T™)  l  (5,  6.67,  673). 


respectively.  Moreover,  the  analysis  and  discussion  of  optimum 
start-up  strategy  is  stated  as  follows: 

•  All  stack  temperatures  are  close  to  the  upper  limit  of  operating 
temperature  as  time  approaches  sufficiently  large,  and  mean¬ 
while  all  stack  voltages  are  kept  at  the  highest  level  without  any 
voltage  drop. 

•  The  stack  start-up  response  can  be  accessed  from  the  zoom-out 
of  each  stack  voltage  response.  At  inlet  jacket  temperatures  T™  = 
723  K,  the  high  steam-to-carbon  ratio  can  reduce  start-up  delay. 

•  The  amount  of  water  flow  strongly  affects  the  impedance  for 
anode  overpotential,  so  the  higher  steam-to-carbon  ratio  would 
induce  the  lower  voltage  output. 

Consequently,  the  system  start-up  as  well  as  the  power  perfor¬ 
mance  can  be  improved  simultaneously  by  the  evaluation  of  inlet 
conditions  (e,  NJ.nxJ.n02 ,  T ™).  Besides,  the  responses  of  some  species 
in  the  DIR-MCFC  with  prescribed  operating  condition  are  shown 
in  Fig.  10(a)-(c).  Notably,  Fig.  10(c)  shows  that  the  high  steam- 
to-carbon  ratio  induces  a  larger  amount  of  hydrogen  but  a  lower 
amount  of  CO.  It  stated  that  the  operating  condition  for  the  stack 
system  would  not  be  prone  to  carbon  monoxide  “poisoning”.  In 
the  presented  preheating  configuration,  it  provides  a  prospective 


alternative  to  improve  the  stack  start-up  operation.  In  fact,  too  short 
warm-up  period  for  a  scale-up  power  system  is  probably  unrealistic 
and  absolute  challenge. 

5.  Conclusions 

In  this  article,  the  new  heat  integration  configuration  of  the 
DIR-MCFC  system  is  proposed.  The  heating  jacket  is  denoted  as 
an  additional  heat  exchanger.  Regarding  the  start-up  operation, 
the  start-up  delay  and  inverse  response  are  investigated  via  the 
dynamic  simulation.  The  manipulation  of  steam-to-carbon  ratio, 
air  flow  rate  and  inlet  jacket  temperature  could  dominate  the  stack 
start-up  operation.  The  steady-state  optimization  algorithm  is  suc¬ 
cessfully  employed.  Under  optimum  air  flow  rate  with  respect  to 
prescribed  inlet  conditions,  the  system  start-up  delay  is  improved 
while  the  maximum  power  output  is  achieved.  Basically,  the  opti¬ 
mum  initial  condition  is  not  unique  and  depends  on  specified 
operating  constraints.  With  aid  of  dynamic  simulation  and  eval¬ 
uation  procedure,  the  system  start-up  delay  is  not  only  reduced 
but  the  satisfactory  power  performance  is  also  guaranteed. 
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Appendix  A. 

Assume  that  both  reversible  reactions  in  Eqs.  (1)  and  (2)  obey 
the  laws  of  chemical  equilibrium,  the  rate  equations  rR  and  rw  are 
formulated  by  [4,18]: 


Ar  exp(-A£R/Rr)[(PCH4PH2o/^25)  -  (P+Pco/Kfe)] 

DEN2 


rw 


Aw  exp  {- AEW / RT)[(PcqPh2o / Ph2  )  -  (Pco2/Kwe)] 
DEN2 


(Al) 

(A2) 


where, 


DEN  =  1  +  BC0  exp (-AHC0/RT)PC0  +  B„2  exp  (- AHh2 /RT)P h2 
+  ficH4  exp  (- AHch4/P1’)Pch4 

+  Bh2o  exp  (- AHH2o/Rr)PH2oPH21  (A3) 


and  equilibrium  constants  are  given  by: 


I<Re  =  1.03  x  1010  exp 


/<We  =  exp 


-3.79762  + 


30.4197- 

4159.54 


27106.2 


(A4) 

(A5) 


Notably,  above  values  of  rate  factors  for  the  internal  reforming 
unit  is  listed  in  Table  1 .  The  rate  for  the  electrochemical  reactions  at 
the  anode  and  cathode  in  Eqs.  (3)  and  (4)  are  given  by  the  Faraday’s 
law  of  electrolysis: 

„  *stackAcelinceii  /Aci 

re  =  2P  (Abj 

Inspired  by  the  dynamic  modeling  of  DIR-MCFC  in  some  litera¬ 
tures  [11,19],  the  mole  balance  equations  for  the  electrochemical 
half  reaction  at  the  anode  are  described  by: 

Xa,i  =  (Na)-1  [N>"  (+,.  -  xaJ)  -  xaJ  (2rR  +  re)  +  R0,i]  (A7) 
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where  the  species  i  eSa  =  {H2,  CH4,  CO,  C02,  H20}  and  Ra  = 
[3 rR  +  rw  -  re,  - rR ,  -  rw,  rw  +  re,  -rR  -  rw] .  The  electrochemi¬ 
cal  half  reaction  at  the  cathode  is  modeled  by: 

*cj  =  (At)-1  [N‘n  (4nj  -  xcJ)  +  1 .5xcjre  +  RCJ]  (A8) 

where  the  species  jeSc  =  {C02,  02},  andRc  =  [-re,  -0.5re].  Since 
the  stack  is  completely  covered  by  heating  jacket,  the  heat  trans¬ 
fer  rate  between  the  stack  and  heating  jacket  is  simplified  by 
[20]: 


Qheat  —  QA 


{Tin_T)_{Tout_T) 


(A9) 


and  the  heating  jacket  is  modeled  by: 


d7T 

dt 


N  (jin  _  'T’OutX  _  ^heat 

"''y  1  hi  ’  mhjCp 


(A10) 


where  T™  and  T°ut  represent  the  inlet  and  outlet  steam  tempera¬ 
ture,  respectively.  The  heat  capacity  of  steam  is  formulated  by: 

Cpj  =  7.256  +  2.298  x  10“3T"ut  +  2.83  x  l(T9(T°ut)2  (All) 

Under  assumptions  (i)-(vii),  the  thermal  model  of  the  stack  is 
described  by  a  first-order  differential  equation: 


dT  _  Ea^g  +  Ec?g  ( T  Tam b)/Rt  ^celdstack  +  Qheat 

dt  msCp  +  {Eah  +  Ech)/T 


(A12) 


where  EaM  =  NiJ'J2ieSaxaA,  Ec,h  =  ^Y.jeScxCA’  E°.g  = 
^£i.sa(  W  -  XaA)  -  and  £c,g  = 

^ZjeSc  ~  X‘Jhi)  ~  Zj,  ScXcfhjJ2RcJ-  The  SPeCifiC 


i  G  Sc 


enthalpies  of  species  at  the  anode  and  cathode,  are  written  by: 


hk  =  h[ef- 


f  (ak  +  bku  +  cku2)  du,  ke  {Sa,Sc}  (A13) 


where  the  value  of  each  hf  and  coefficients  for  each  component 
are  listed  in  Table  2.  Finally,  the  performance  of  fuel  cell  is  eval¬ 
uated  by  the  cell  voltage  under  a  specified  current  demand.  The 
cell  voltage  of  DIR-MCFC  consists  of  equilibrium  cell  potential,  cell 
overpotentials  and  cell  resistance: 


VceH  —  Feq  i stack  (^?ohm  +  V a  +  Tie) 


where  the  equilibrium  cell  potential  is  described  by: 


„  2.44  xlO5 

Feq  —  - 


(472.45  +  0.04T)  T  RT  Pa,H2Pc,co2P°o2 
2  F  +  2  F  Pa,  C02  Pa,  h2o 

(A15) 


Cell  overpotentials  are  based  on  the  following  empirical  equa¬ 
tions  [21]  to  describe  anode  and  cathode  overpotentials: 


Tja  =  2.27  x  10  9  exp 


/  6435 


rjc  =  7.505  x  10“10  exp 


n-0.42  n-0.17  n-1 
1  a,H2  a,C02  a,H20 


,-0.43  p-0.09 

02  1  c,C02 


(A16) 


Cell  resistance  which  contributes  to  the  ohmic  loss  is  expressed 
with  the  Arrhenius  equation  [22]: 


^7ohm  —  9.5  x  10  exp 
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